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The reduced-mass dependence of relativistic and radiative effects in simple muonic bound sys- 
tems is investigated. The spin-dependent nuclear recoil correction of order (Za)'^ fi^ /m% is evaluated 
for muonic hydrogen and deuterium, and muonic helium ions (/i is the reduced mass and ttijv is 
the nuclear mass). Relativistic corrections to vacuum polarization of order a(Za)'^fi are calculated, 
with a full account of the reduced-mass dependence. The results shift theoretical predictions. The 
radiative-recoil correction to vacuum polarization of order a{Za)^ \n^{Za)fj,^ /mj^ is obtained in 
leading logarithmic approximation. The results emphasize the need for a unified treatment of rel- 
ativistic corrections to vacuum polarization in muonic hydrogen, muonic deuterium and muonic 
helium ions, where the mass ratio of the orbiting particle to the nuclear mass is larger than the 
fine-structure constant. 

PACS numbers: 12.20.Ds, 36.10.Ee, 14.20.Dh, 31.30.jf, 31.30.jr 
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I. INTRODUCTION 

In muonic hydrogen and muonic deuterium, the mass 
ratio ^AT = m^/mN of the orbiting particle (muon mass 
771^) to the mass of the atomic nucleus rriN is not really 
small against unity. It evaluates to 



^p^—^ =0.112609. 



Cd = — ^ ==0.0563327. 
md 
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where the latest recommended values of the masses have 
been used [l|. For muonic helium ions, we have 



^Ho 



mHcS 

771, 



:= 0.0376223 ...« ^ , (Ic) 



^Hoi = — ^ ^ 0.0283465 . . . « :^ . (Id) 

In all cases, S,n is larger than the fine-structure constant 
a « 1/137.036 that governs the relativistic and quan- 
tum electrodynamic (QED) effects. Consequently, the 
reduced-mass dependence of all quantum electrodynamic 
(QED) effects that influence the spectrum must be taken 
into account exactly, i.e., to all orders. In calculations, 
one must flrst take into account ^n (if possible) to all 
orders, before advancing to the next order in the Za- 
expansion; otherwise the higher-order effects in Za will 
be shadowed by the unknown reduced-nrass dependence 
of lower-order terms in the Za-expansion. 

Hence, particular emphasis has been laid in Ref. [2| on 
the correct treatment of the reduced-mass dependence of 
all relativistic and QED corrections. The statement made 
in the text preceding Eq. (17) in Ref. |2[, which says "the 
external field approximation docs not give an accurate re- 
sult," can hardly be overemphasized. Here, the external 
field approximation refers to the Dirac equation, which is 
appropriate for heavy muonic atoms where the parameter 
Za (with Z denoting the nuclear charge number) is much 
larger than the mass ratio m^/r?7,7v, where mjv is the mass 



of the heavy nucleus. Even a tiny conceivable error in the 
handling of, say, the reduced- mass dependence of the one- 
loop vacuum polarization (VP) shift in muonic hydrogen 
could drastically influence the comparison of theory and 
experiment: the current discrepancy [3| of theory and ex- 
periment for the muonic hydrogen Lamb shift amounts to 
roughly 0.3 meV, which is about one part per thousand 
of the leading vacuum-polarization contribution and thus 
smaller than a conceivable additional reduced-mass cor- 
rection to the leading VP effect of relative order ^^. 

In comparison to previous studies on heavy muonic 
atonrs and ions (excellent theoretical overviews are pro- 
vided in Refs. 0, Q), the magnitude of the mass ratio 
is the main characteristic property of muonic hydrogen 
and deuterium. In this article, we thus revisit the pre- 
cise treatment of the vacuum-polarization contribution 
to the Lamb shift in muonic hydrogen (fiil) and muonic 
deuterium (/^D), as well as muonic helium ions (/xHe"^ and 
/iHe'^), with full account of the two-body structure of the 
bound system. Starting from the nonrelativistic Hamilto- 
nian (Sec.|ll]), we proceed to discuss the nuclear-spin de- 
pendent terms in the Breit Hamiltonian (Sec. lIII|) . before 
proceeding to the radiatively corrected Breit Hamiltonian 
(Sec. IIV[) and the radiative- recoil correction (Sec. FV)) . 
Conclusions are drawn in Sec. I VII 



II. NONRELATIVISTIC HAMILTONIAN 

The nonrelativistic /iH Hamiltonian is separable, and 
the nonrelativistic (Schrodinger) Hamiltonian in the 
center-of-mass system, where the muon and the nuclear 
particle carry opposite momenta p and —p, respectively, 
reads (in natural units, h = c = eo = 1), 
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This equation can be solved exactly in terms of 
Schrodinger eigenstates. The nonrelativistic spinor wave 



functions for the 25*1/2 and 2Pi/2 states are exact eigcn- 
states of H and read, explicitly, 

^25(r) = ^-^^^ (2 - Za^^r)c~'^^-^^ x'A{r) , (3) 



V'2Pi/2 (^) 
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where M = ±i is the magnetic projection, xi^('LLis the 
standard two-component spin-angular function [6|, and 
><• = (—1)^+^+1/2 is the Dirac angular quantum number. 
The reduced-mass dependence of the wave functions in 
Eq. ([3]) is exact. 

The one-loop vacuum-polarization potential V^p can 
be expressed in terms of the action of a linear operator 
iiT on a screened Coulomb potential Uvp as follows, 

V^p{r) = K[vvp{mep;r)] , Uvp(A;r) = e"^*" , (5) 

r 

with 



out the integration over the spectral parameter p analyt- 
ically, with the result 

(2Pi/2|Kp|2Pi/2> - (25i/2|Kp|25i/2> 
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for the Lamb shift difference of the leading VP en- 
ergy correction. For reference, the 2Pi/2~2S'i/2 differ- 
ence of the leading nonrelativistic vacuum polarization 
effect is 205.0073 mcV for ^H, 227.6346 mcV for /iD, 
1641.885 meV for ^He^, and 1665.772 meV for ^He^. 
The latter value differs by 0.010 meV from the value of 
1665.782 meV given in Eq. (10) of Ref. [il|; the differ- 
ence probably is due to updated physical constants used 
in our calculation (see also Ref. [l|). 
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III. BREIT HAMILTONIAN AND 
BARKER-GLOVER TERMS 



where rrie is the electron mass. In the following, we often 
use the identification A = m.e p and define the ratio 



A 



N 



{Zap) 



(7) 



which evaluates to /3p = 0.7373836 ... for /xH and Pd = 
0.7000861 ... for ^D (proton and deuteron nuclei, respec- 
tively). For muonic helium ions, the values are /^hc-^ = 
0.3438429 ... and f3^^i = 0.3407691 .... (We here refrain 
from assigning a subscript to the reduced mass p, even 
though it of course depends on the nucleus N, because 
the symbol pN is reserved, canonically, for the nuclear 
magnetic moment.) We then use the exact nonrelativis- 
tic unperturbed wave functions defined in Eq. Q and 
calculate the leading VP energy shifts as 



(25*1/2 1 Kp 1 251/2) 
and 
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A numerical evaluation of these compact expressions 
is found to be in agreement with the literature (see 
Refs. l^iSSQ) and confirms that the reduced-mass de- 
pendence of the leading VP effect is correctly described 
by Schrodinger wave functions scaled with the reduced 
mass of the system. It is even possible [3, ll^ to carry 



The Breit Equation and the corresponding Hamilto- 
nian follow from the Bethe-Salpeter equation in the limit 
on an instantaneous interaction kernel |12| and describe 
the bound states of general two-body systems of arbitrary 
mass ratio ^at, including higher-order relativistic correc- 
tions [i3|. For the 2Pi/2-25i/2 Lamb shift in muonic 
bound systems, the relevant terms in the Brcit Hamilto- 
nian read ((5/ = 1 for half-integer and 5/ = for integer 
nuclear spin, see |14| ') 
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where we use the summation convention for the super- 
scripts i and j which denote the Cartesian components 
of the position and momentum operators. Using the re- 
lations 



VM - ) =-47r(5^(r) 



and 



V'V^ 



-f47rJ^(r), 



(11) 



(12) 



one may transform SH^ to a more symmetric form, 
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After some algebra, the expectation values of the eigen- 
states given in Eq. ^ of the Breit Hamiltonian read 



/OQ IXTJlOQ \ (V \4 5 + ^Ar(ll + 13^Ar) 

\25i/2 \0il\ 2bi/2) = - [Za) fi- 
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and the 2Pi/2-2S'i/2 difference (2Pi/2 is energetically 
higher) amounts to 



L(2Pi/2-25i/2) 
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The Barker-Glover [15[ correction L given in Eq. ([15]) 
evaluates to 0.05747 meV for /xH and to 0.12654 meV for 
/iHe'^, in full agreement with the literature [Eq. (46) of 
Ref. Q]. Because the zitterbewegung term is absent for 
the spin-1 deuteron nucleus [Tj] and for the spin-zero al- 
pha particle, the shift evaluates to L = 0.06722 meV for 
^D and to L = 0.29518 meV for fiUe'^ [d. Eq. (61) of 
Ref. [HI and Eq. (10) of Ref. [H]]. It constitutes a nu- 
clear spin-dependent recoil correction to the Lamb shift, 
which is essential for the correct description of the niuonic 
isotope shift. Equation p3|) is exact to all orders in ^jy. 



IV. RADIATIVELY CORRECTED BREIT 
HAMILTONIAN 



TABLE I. Detailed breakdown of the first-order and second- 
order individual contributions SE^ and SE) to the rela- 
tivistic Breit correction SE^p of vacuum polarization for fiYL, 
fiD, and muonic helium ions. All units are meV. 



1^ 
SEi'^ 

SEi'^ 



SEW 






5£(2 



5£vp 



7^ 
SEi'^ 
SEi'' 
SEi'^ 



5£(i) 



Ie^ 

SEi'^ 

se\'^ 

5Ef^ 



/iH 



/xD 



^He'^' 



AtHe* 



2Pi 



/2 



eV] 



-0.000558 

0.000064 

-0.000290 

-0.002026 



-0.000679 

0.000038 

-0.000181 

-0.002303 



-0.020331 

0.000467 

-0.004587 

-0.085970 



-0.020970 

0.000360 

-0.003584 

-0.087587 



-0.002811 -0.003125 -0.110421 -0.111781 



-0.001124 
0.0 

-0.000269 
-0.001283 



-0.001545 
0.0 

-0.000177 
-0.001521 



-0.099980 
0.0 

-0.008427 
-0.093497 



-0.105132 
0.0 

-0.006624 
-0.095762 



-0.002676 -0.003243 -0.201904 -0.207518 



-0.005486 -0.006368 -0.312324 -0.319300 



2Si 



/2 



eV] 



0.029112 
-0.001928 
-0.002280 

0.0 



0.034636 
-0.001142 
-0.001416 

0.0 



0.846700 
-0.014512 
-0.032734 

0.0 



0.872150 
-0.011243 
-0.025535 

0.0 



0.024904 0.032078 0.799454 0.835372 



-0.084996 
0.044911 

-0.009064 
0.0 



-0.108282 
0.053594 

-0.005539 
0.0 



-2.875794 
1.361115 

-0.106444 
0.0 



-2.995690 
1.402803 

-0.082889 
0.0 



5£(2) 


-0.049149 -0.060227 -1.621122 -1.675776 


5£vp 


-0.024245 -0.028149 -0.821668 -0.840404 


2^1/2-251/2 [meV] and comparison to other work 



ASvp (this work) 0.018759 0.021781 0.509344 0.521104 

(Ref. [17]) 0.0169 

(Ref. [nj) -0.202 

(Ref. [lel)" 0.0169 0.0214 0.495 0.508 

° A conceptually different approach is used in Ref. [16|]. 



The massive Breit interaction uses a strictly static 
timelike photon propagator component 



transforms, the most interesting of which is related to the 
interaction [cL Eq. (83.13) of jl| and Sec. 2 of RcL O], 



and spatial components 



G^M = - 
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The spatial components are no longer transverse. One 
then follows the standard derivation of the Breit interac- 
tion given in Chap. 83 of Ref. [T^ but has to avoid pit- 
falls. The derivation necessitates the evaluation of Fourier 
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For A = 0, the Fourier transforin of this expression with 
respect to (f gives the term 5H^ in Eq. ([TU| . For a massive 
photon, we find 



d^g 

(27r)3 



Uip.q^Xy"!-^ ^ 5v2{r) + 6v^{r), 



(19) 



where dv2{r) and Sv^^r) contribute to the Brcit potential 
Svvp for massive photon exchange, 

(5wvp = K[5vi + 5v2 + Sv3 + Sva] , (20) 

where 5vi(r) depends on the nuclear spin, 

Za ( 1 
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and the niomentum operators act on the ket state in 
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whereas the spin-orbit coupling is modified to 



6v4 = Za 



( 1 



1 



\47Ti^ 2vi^^mN 
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(1 + \r) 



a-L. (23) 



In the terms dv2 and Sv^, all the momentum operators 
act on the "incoming" wave function (Dirac ket state), 
and the Hamiltonian may be used for the evaluation of 
diagonal matrix elements. For off-diagonal elements, it is 
helpful to symmetrize Sv2 and Sv^ with respect to outgo- 
ing and incoming momenta, effectively replacing terms 
of the form /(r)ir • p by the commutator i[/(r)f*, p] 
and terms of the form J*-' (r) p* fp by the anticommuta- 
tor i {p-'{r), P^ P'}- In a second step, using the relation 
i {A^,B} = AB A + \[A, [A,B]], one obtains an even 
more symmetric form, with 



5wi = 5vi , 5wi = (5w4 



5W2 = - 



5w?^ 



ZaVe 
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AmnTriNr 



Xr\ 
1-^ , 
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(24b) 



p [ — ^ ^ — rV-' p' , (24c) 
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Svvp = K[Swi + Sw2 + Sw3 + Sw4 



(24d) 



The terms S'W2 and 6W3 are used in Eq. (21) of Ref. [3- 
The a{Za)^ 11 relativistic reduced-mass correction to vac- 
uum polarization then is the sum of four first-order per- 
turbations 5E\ and four second-order terms SE- , 



SE^p = <5i?(i) + <5i?(2) = ^ 5ij(i) + ^ SEf , (25a) 

SeI'^ ^K[{nlj\dw^\n£j)] , 
6Ef ^2K[{ne,\6H,\S^„i^)] , 



(25b) 
(25c) 



where \Sipne) is the wave function correction due to VP, 

1 



\Hne,) 



Eril — H 



Wvp \nlj) 



(26) 



Using a generalization of techniques outlined in Ref. [20| , 
the perturbation Sipntj can be evaluated analytically. The 
detailed expressions for the reduced Green functions (in- 
dicated by a prime) of the 2Si/2 and 2Pi/2 states have 
been given in Eqs. (23) and (24) of Ref. 0. All individual 
contributions are listed in Table HI in order to facilitate 
a numerical comparison with independent calculations. 
For jiH, we obtain a result of AE^p — 6E^p{2Pij2) — 
6E^p{2Sii2) = 0.018759 meV. This result is not in per- 
fect agreement with published values [2, [la, [l3|- For 
comparison, the result indicated in Eq. (25) of Ref. [2| 
reads 0.059 meV; and in Eq. (25) of Ref. [13 a result of 
0.0169 meV has been indicated. In Table 1 on page 8 of 
Ref. [16[, a numerically equivalent result of 0.0169 meV 
is given. (We note that Ref. [16] contains many unnum- 



bered tables; the referenced table is numbered). The ma- 
trix elements of the relativistic recoil operator given in 
Eq. (7) of Ref. [IQl are evaluated using unperturbed wave 
functions. All values given in TableUare nonperturbative 
in the mass ratio and take the wave function correction 
into account. A precise comparison of individual contri- 
butions to the approach of Ref. [ia| is not possible at 
present. As evident from Table HI there are quite signifi- 
cant differences with published values for /zHe^: e.g., the 
entries in Eqs. (26)-(29) and Eq. (41) of Ref. [Ill add up 
to a correction of — 0.202 meV for the 2Pi/2^2S'i/2 Lamb 
shift in /iHe^, whereas we obtain -1-0.521 meV. 

A very important question concerns the verifiability 
of the results. In self-energy calculations [21[, a cross- 
check of the calculation consists in the cancellation of an 
overlapping parameter that separates different momen- 
tum and energy regions of the physical process. For VP 
effects in muonic systems, no such checks are immedi- 
ately available. Here, we note that the entries for the 
first-order matrix elements in Table U for /iHe'* are in 
full agreement with the results given in Eqs. (26)-(29) 
of Ref. [ll|. For the matrix elements needed for 6E^^\ 
the limit as A — >■ of the matrix elements (n£j \5'Wi\ n£j) 
can be verified independently, and the calculation can 
otherwise be performed analytically, with ease. For the 
matrix elements needed in the evaluation of the second- 
order effects SE^^\ we can verify the first few terms in 
the asymptotic limit as A — )■ 0, using the relation 



2 ( n£. 



SH 



1 



= 2 ( nl. 
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n£j ) X^+0{X^)^ 



In deriving this relation, the Hellmann-Feynman theo- 
rem is useful for the zeroth-order term in A. The wave 
function perturbation in the term of order A^ can be eval- 
uated analytically. 





FIG. 1. Feynman diagrams for the radiative-recoil correction 
in two-body muonic bound systems. The three given diagrams 
correspond to the vacuum-polarization insertion in the seag- 
ull and two-photon exchange and lead to the leading double 
logarithm given in Eq. 



V. RECOIL CORRECTIONS TO VACUUM 
POLARIZATION 

Beyond the radiative modifications of the static Breit 
Hamiltonian, the recoil correction to vacuum polarization 
can be obtained by the insertion of vacuum polarization 
loops into the Salpeter recoil correction |23 - |24| . The re- 
coil correction is the sum of four terms [2J1 ; two of these 
(low- and middle-energy part) describe the frequency- 
dependent part of the Breit interaction, beyond the static 
Breit Hamiltonian given in Eq. (jlO[) . and two further 
terms (seagull and high-energy part) correspond to two- 
photon exchange. 

The seagull term corresponding to Fig. [1] (left), with 
a vacuum polarization insertion in the exchange photon, 
leads to the integral 



5Es^ 



5'^ - 




(28) 



nl,Wi^'+^-')-'^\nij 



where Wi = y k^ + X^ is the frequency of the massive 
photon in the vacuum-polarization loop. An ultraviolet 
cutoff A is introduced via multiplication of the integrand 
by a multiplicative regularization factor -=5 



The auxiliary parameter A cancels when the high-energy 
part from two-photon exchange [see Figs. [TJ middle, 
and Fig.IU right] is added to the result (see also Ref. j24|). 
From the integral (^5]) . we extract a leading double log- 
arithmic correction, 



6Es 



Aa{ZaffiHe, 



■ In^ {4:Za0 



(29) 



which is nonvanishing only for S states {£ — 0). This 
correction evaluates to 0.0003 meV for the 2Pi/2-25i/2 
Lamb shift in ^H, 0.0002 meV for ^D, 0.0072 meV for 
fiHe^, and 0.0056 meV for ^He''^. Because subleading log- 
arithmic terms, and nonlogarithmic terms are missing, 
the theoretical uncertainty of the results in Eq. ([29)) 
should be taken as 100 % of the leading logarithmic cor- 
rection calculated here. 



VI. CONCLUSIONS 

Our theoretical investigations are motivated by the ne- 
cessity to shed light on the recently observed discrepancy 
of theory and experiment in /xH (see Ref. [3| . By an ex- 
plicit evaluation of the matrix elements of the two-body 
Breit Hamiltonian, we obtain the nuclear-spin dependent 
recoil contributions to the Lamb shift in fjH and /nD given 
in Eq. ([T5|). and confirm that the results are exact in 
the mass ratio, so that the existence of further recoil 
corrections [l5| can be ruled out at order (Za)^. The 
calculation of the relativistic reduced-mass corrections 
to vacuum polarization using the massive Breit Hamil- 
tonian is shown to involve a nontrivial nuclear-spin de- 
pendent term [see Eq. (PT|) ]. Our detailed numerical in- 
vestigation (see Table HI) slightly decreases the observed 
experimental-theoretical discrepancy [3| (in contrast to 
a recent investigation [25|, where the authors obtain an 
increase of the discrepancy, based on a treatment which 
is perturbative in the mass ratio). 

A detailed breakdown of the relativistic corrections to 
vacuum polarization, including the reduced-mass correc- 
tions is given in Table |T] for muonic hydrogen, muonic 
deuterium and muonic helium ions. For muonic hydro- 
gen, the sum of the entries in rows 3 and 19 of the theory 
in the supplemental material of Ref. (3| , minus the entry 
in row 1 of the same supplemental material, amounts to 
(205.0282 - 0.0041 - 205.0074) meV = 0.0167 meV; this 
is close to the result indicated in Table 1 of Ref. [16| . 
which is 0.0169 meV. In Refs. 0, 0, [la| , the second entry 
in the mentioned combination (—0.0041 me V) has been 
referred to as a "recoil correction to vacuum polariza- 
tion" , whereas we here refer to the effect as a relativistic 
correction to vacuum polarization with a proper account 
of the reduced-mass dependence. Our approach is non- 
perturbative in the mass ratio ^n and isolates the terms 
of order a (Za)'^, while treating the two-body aspects of 
the problem to all orders. 

Our calculations lead to significant shifts of theoretical 
predictions for /xHe'^ with respect to published values (ex- 
periments are planned for the near future). Specifically, 
for /iHe"*, our nuclear-spin dependent Barker-Glover type 
correction L of 0.295 meV differs from the value of 
0.074 meV given in Refs. [H Ull by -H0.221meV. For 
fiHe'^ , our result for the relativistic correction to vacuum 
polarization, with a full account of the reduced-mass de- 
pendence, reads as 0.521 meV for the 2P-2S difference, 
to be compared with a value of — 0.202 meV given in 
Ref. [Il|. This leads to a total upward shift of theoret- 
ical predictions for the 2Pi/2-2Sif2 Lamb shift in /iHe^ 
by [0.521 - (-0.202) + 0.221 - 0.010]meV = +0.934 meV 
relative to Ref. [ll| (where we add the small correction of 
the reference value of the leading VP correction) and by 
(+0.221 + 0.013) meV = 0.234 meV relative to Ref. [3 
(where we add the difference in the relativistic correction 
to vacuum polarization from Table H]). 

The radiative-recoil correction obtained in Eq. (|29p is 
numerically small; however, this two-loop bound-state 



correction has traditionally been one of the most elu- 
sive effects in bound-state quantum electrodynamics for 
two-body systems. Its calculation in leading logarithmic 
approximation helps to determine the overall uncertainty 
of theoretical predictions with regard to the conceptually 
involved higher-order recoil corrections to VP, given by 



the two-body nature of the bound system. 
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